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The evaluation of the hadronic contribution to the muon magnetic anomaly is reviewed, including 
a new estimate using precise results on the tt + tt~ spectral function from the KLOE Collaboration. 
It is found that the KLOE data confirm to some extent the previous e + e — annihilation data in 
this channel, and accentuate the disagreement with the isospin-breaking-corrected spectral function 
from r~ — > Tr~n°v T decays. Correcting for the empirical difference in the mass of the charged and 
the neutral p locally improves, but does not resolve this discrepancy. A preliminary reevaluation 
(including the KLOE data) of the e + e — -based Standard Model prediction of results in a deviation 
of 2.7 standard deviations from the BNL measurement. 



1 Introduction 

Hadronic vacuum polarization (HVP) in the 
photon propagator plays an important role in 
many precision tests of the Standard Model. 
This is the case for the muon anomalous mag- 
netic moment a M = (g^ — 2)/2, where the 
HVP component is the leading contributor to 
the uncertainty of the Standard Model pre- 
diction. The HVP contribution is computed 
by means of a dispersion relation as an inte- 
gral over experimentally determined spectral 
functions. It is the property of this disper- 
sion relation that the 7T7t spectral function 
provides the major part of the total HVP 
contribution, so that the experimental effort 
focuses on this channel. 

Spectral functions are directly obtained 
from the cross sections of e + e~ annihilation 
into hadrons. The accuracy of the calcu- 
lations has therefore followed the progress 
in the quality of the corresponding dateH 
Because the latter were not always suit- 
able, it was deemed necessary to resort to 
other sources of information. One such pos- 
sibility was the use of the vector spectral 
functions^ derived from the study of hadronic 
r decayJ^for the energy range less than m T ~ 
1.8GeV/c 2 . For this purpose, the isospin ro- 
tation that leads from the charged r to the 
neutral e + e~ final state has to be thoroughly 



corrected for isospin-breaking effects. 

Also, it was demonstrated that essen- 
tially perturbative QCD could be a pplied to 
energy scales as low as thus of- 

fering a way to replace poor e + e~ data in 
some energy regions by a reliable and precise 
theoretical prescriptiorl ^ l ^ l ^ l ^^^^^ l 

Detailed reanalyses including all avail- 
able experimental data have been published 
in RefsP^IE] (s ee also the preliminary re- 
sults given in RefsHnnsi), taking advantage 

of precise results i n th e 7T7t channel from the 
CMD-2 experiment and from the ALEPH 
analysis of r decays^! and benefiting from a 
more com plete tr eatment of isospin-breaking 
correction J ^ 1 ^ 1 It was found that the e + e~ 



and the isospin-breaking-corrected r spectral 
functions do not agree within their respec- 
tive uncertainties, thus leading to inconsis- 
tent predictions for the lowest-order hadronic 
contribution to a^. The dominant contri- 
bution to the discrepancy stems from the 
7T7T channel with a difference of (—11.9 ± 
6.4 exp ±2.4 rad ±2.6 su(2) (±7.3 to tai)) x KT 10 , 
and a more significant energy-dependent de- 
viation. When compared to the world av- 



a The problem between t and e + e — data is more not- 
icable when comparing the r — ► nir°u branching frac- 
tion with the prediction obtained from integrating the 
corresponding isospin-breaking-corrected e+e — spec- 
tral function. Here, the function under the integrand 
is less selective than it is the case for the HVP contri- 
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erage of the muon magnetic anomaly mea- 
surements, dominated by the results from the 
BNL experiment!^, 



(11 659 208.0 ±5.8) x 10" 



(1) 



the respective e + e~ and r-based predictions 
disagreed at the level of 2.5 and 1.3 stan- 
dard deviations, when adding experimental 
and theoretical errors in quadrature. 

This summer, new data on the mr spec- 
tral function in the mass region between 0.60 
and 0.97 GeV/c 2 were presented by the KLOE 
Collaboration , using the — for the purpose 
of precision measurements — inn ovative tech- 
nique of the radiative returrP^I. The statisti- 
cal precision of these data by far outperforms 
the Novosibirsk sample, but the systematic 
errors are about twice as large as those ob- 
tained by CMD-2. New data using the same 
technique have been published by the BABAR 
Collaboration^! on the ir + TT~ir° final state. 
They unveil a larger cross sections and a res- 
onant peak at around 1.6 GeV/ c 2 that was 
missed by the previous DM2 measurement!*^! 
The BABAR data are not (yet) used in the 
preliminary reevaluation of the lowest-order 
HVP contribution given here. 6 

2 Muon Magnetic Anomaly 

It is convenient to separate the Standard 
Model (SM) prediction for the anomalous 
magnetic moment of the muon into its dif- 
ferent contributions, 
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and where ag ED = (11658 472.0 ± 0.2) x 
is t he pur e electromagnetic contri- 
and references therein), 
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bution to a u , leading to a discrepancy of 2.9 standard 
deviation*"*! 

^The correction to aJ^ ad ' L when using the BABAR. 
data for this mode is of the order of +1 X 10 — 10 . 



a had,LO j g foe lowest-order HVP contribu- 
tion, a £ ad ' HO = (_ 10 o ± fj 6) x 10~ 10 is 
the corresponding higher-order parlPEl and 
a wcak = (15.4 ±o.l ±0.2) x 10" 10 , where the 
first error is the hadronic uncertainty and the 
second is due to the Higgs mass range, ac- 
counts for corrections due to exchange of the 
weakly interacting bosons up to two loops^l 

For the light-by-light (LBL) scattering part, 
fl had,LBL 5 we uge the value ( 12 o ± 3.5) x 

10~ 10 from the latest evaluation^!, slightly 
corrected for the missing contribution from 
(mainly) the pion box. 

Owing to unitarity and to the analytic- 
ity of the vacuum-polarization function, the 
lowest order HVP contribution to a., can be 
computed via the dispersion integral"^ 



liad.LO 
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where K(s) is a well-known QED kernel, and 
(s) denotes the ratio of the "bare" cross 
section for e + e~ annihilation into hadrons to 
the pointlike muon-pair cross section. The 
bare cross section is defined as the mea- 
sured cross section corrected for initial-state 
radiation, electron-vertex loop contributions 
and vacuum-polarization effects in the pho- 
ton propagator. However, photon radiation 
in the final state is included in the bare cross 
section defined here. The reason for using the 
bare (i.e., lowest order) cross section is that 
a full treatment of higher orders is anyhow 
needed at the level of a^, so that the use of 
the "dressed" cross section would entail the 
risk of double-counting some of the higher- 
order contributions. 

The function K(s) ~ 1/s in Eq. (0J gives 
a strong weight to the low-energy part of the 
integral. About 91% of the total contribu- 
tion to a had . LO is accumulated at center-of- 
mass energies y/s below 1.8 GeV and 73% of 
a had,LO j s coverec ] by the tttt final state, which 
is dominated by the p(770) resonance. 
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Figure 1. Relative comparison of the vt+tt spectral functions from e + e -annihilation data and isospin- 
breaking-corrected r data, expressed as a ratio to t he t spect ral function. The shaded band indicates the 
errors of the r data. The e+e~ data are from KLOrBH CMD-JIZJ CMD, OLYA and DM1 (references given 
m The right hand plot emphasizes the region of the p peak. 



3 The Input Data 

A detailed compilation of all the experimen- 
tal data used in the evaluation of t he dis- 
persion integral J3J is provided in 
Also discussed therein is the corrective treat- 
ment of radiative effects applied to some of 
the measurements. The r spectral function 
is obtained by averaging the re sults from 
ALEPlP, CLECJSS and OPAlESI which ex- 
hibit satisfactory mutual agreement. 

A comparison of the e + e~ — > tt + it~ data 
and the corresponding r spectral function, 
represented as a point-by-point ratio to the 
r spectral function is given in Fig.^ Several 
observations can be made. 

• A significant discrepancy, mainly above 
the p peak is found between r and the 
e + e~ data from CMD-2 as well as older 
data from OLYA. 

• Overall, the KLOE data confirm the 
trend exhibited by the other e + e~ data. 

• Some disagreement between KLOE and 
CMD-2 occurs on the low mass side 
(KLOE data are large), on the p peak 
(KLOE below CMD-2) as well as on the 
high mass side (KLOE data are low). 



At this stage, the r spectral function has 
not been corrected for a p ossible p~ —p° mass 
and width splitting In contrast to ear- 

lier experimental^ and theoretical resultJ^, 
a combined pion form factor frPO to the new 
precise data on e + e _ and r spectral functions 
leads to m p - - m p0 = (2.3 ± 0.8)MeV/c 2 , 
while no significant width splitting c is ob- 
served within the fit error of 1.7MeV/c 2 . 

Considering the mass splitting in the 
isospin-breaking correction of the r spectral 
function tends to locally improve (though not 
restore) the agreement between r and CMD- 
2 data, leaving an overall normalization dis- 
crepancy. Increasing the r„- — r p o width 
splitting by +3 MeV/c 2 improves the agree- 
ment between r and KLOE data in the peak 
region, while it cannot correct the discrepan- 
cies in the tails. Note that a correction of 
the mass splitting alone would increase the 
discrepancy between the r and e + e~-based 



c Note that if the mass difference is to be taken 
as an experimental fact, a larger width differ- 
ence would be expected . U sing a chiral model 
of the p resonanc one has To = 

r p _(m p0 /m p -) 3 (/? /7?_) 3 + Ar EM , where Ar EM 
is the width difference from electromagnetic decays. 
This leads to a total width difference of (2.1 ± 
0.5) MeV/c 2 tha t is marginally consistent with the ob- 
served valuJ2£( 
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Figure 2. Comp ariso n of the result JgJ with the 
BNL me asur ement^ Also given are our previous 
estimateJi^J where the triangle with the dotted er- 
ror bar indicates theTj-based result, as well as the 
estimates from Refs^^E^^SJ not yet including the 
KLOE data. 



results for a^ ad,L °. 



,had,LO 



During the previous evaluations of 
, the results using respectively the r 
and e + e~ data were quoted individually, but 
on the same footing since the e + e~-based 
evaluation was dominated by the data from 
a single experiment (CMD-2). The confirma- 
tion of this discrepancy by KLOE discredits 
the r-based result for the use in the disper- 
sion integral until a better understanding of 
the dynamical origin of the observed effect 
is achieved. This is a challenging problem, 
which may itself turn out to be of fundamen- 
tal importance. 

4 Results 

The inclusion of the KLOE tttt data decreases 
the contribution from this mode froirM^ 1 
(450.2±4.9±1.6 rad ) x 1CT 10 to (448.3 ± 4.1 ± 
l-6 ra d) x 10 -10 for the energy interval between 
0.5 and 1.8 GeV. Note that the additional 
systematic error due to radiative effects orig- 
inates from the energy regions not covered 
by the recent KLOE and CMD-2 measure- 
ments, where a full treatment of radiative 



corrections is applied. The preliminary esti- 
mate of the integral J3J given below includes 
one additional improvement with respect to 
Ref.El perturbative QCD is used instead of 
experimental data in the region between 1.8 
and 3.7 GeV, where non-perturbative contri- 
butions to integrals over differently weighed 
spectral functions were found to be smalH 
This results in a reduction of a c ^ d ' lj0 by 
— 1 x 10~ 10 . All other contributions to the 
disper sion integral are equal to those defined 
in RefPD. 

The e + e~ -based result for the lowest or- 
der hadronic contribution is 



a had,LO = (693.4±5.3±3.5 rad )xl(r 



(5) 



where the second error is due to our treat- 
ment of (potentially) missing radiative cor- 
rections in the older datcM^I Adding to 
this the QED, higher-order hadronic, light- 
by-light scattering, and weak contributions 
given in Section [21 one finds 



„SM 



(11 659 182.8 ±6.3 h ad,LO+HO 

± 3.5had,LBL ± 0.3qed+Ew) x 10~ 



(6) 



This value can be compared to the 
present measurement adding all errors 
in quadrature, the difference between exper- 
iment and theory is 



cxp 



,SM 



(25.2 ±9.2) x 10 



-10 



(7) 



which corresponds to 2.7 "standard devia- 
tions" (to be interpreted with care due to 
the dominance of systematic errors in the SM 
prediction). A graphical comparison of the 
result © with previous evaluations and the 
experimental value is given in Fig. [21 

5 Conclusion and Perspective 

In spite of the new and precise data 
on the two-pion spectral function from 
the KLOE Collaboration, the lowest or- 
der hadronic vacuum-polarization contribu- 
tion remains the most critical component in 
the Standard Model prediction of a„. The 
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central piece of information provided by the 
present KLOE data (a reanalysis with better 
quality data and a refin ed s tudy of systematic 
effects can be expected^) is that they con- 
firm the discrepancy between the r data and 
e + e~ annihilation observed in this channe^^. 
This said, we point out that there also oc- 
curs disagreement between KLOE and CMD- 
2 data in some of the energy regions. 

An empirical isospin-breaking correction 
of the p resonance lineshape (mass and 
width) improves but does not restore the 
agreement between the two data sets. It is 
a consequence of this confirmation that, un- 
til the CVC puzzle is solved, only e + e~ data 
should be used for the evaluation of the dis- 
persion integral. Doing so, and including the 
KLOE data, we find that the Standard Model 
prediction of a M differs from the experimental 
value by 2.7 standard deviations. 

We are looking forward to the forthcom- 
ing results on the low- and high-energy two- 
pion spectral function from the CMD-2 Col- 
laboration. These data will help to signifi- 
cantly reduce the systematic uncertainty due 
to the corrective treatment of radiative ef- 
fects, often omitted by part by the previous 
experiments. 

The initial-state-radiation program of 
the BABAR collaboration has already proved 
its performance by publishing the spectral 
function for 7T + 7r~7r° (and soon for 2tt + 2tt^), 
while results for the two-pion final state are 
expected. 



Acknowledgements 

I am indebted to Arkady Vainshtein and 
Graziano Vcnanzoni for informative discus- 
sions, and to Bill Marciano for pointing out a 
mistake in the citation of the QED contribu- 
tion. I thank Michel Davier, Simon Eidclman 
and Zhiqing Zhang for the fruitful collabora- 
tion. 



References 

1. S. Eidelman and F. Jegerlehner, Z. Phys. 
C67, 585 (1995). 

2. R. Alemany, M. Davier and A. Hocker, 
Eur. Phys. J. C2, 123 (1998). 

3. ALEPH Collaboration (R. Barate et al), 
Z. Phys. C76, 15 (1997). 

4. ALEPH Collaboration (R. Barate et al), 
Eur. J. Phys. C4, 409 (1998). 

5. OPAL Collaboration (G. Abbiendi et 
al), Eur. Phys. J. C7, 571 (1999). 

6. A.D. Martin and D. Zeppenfeld, Phys. 
Lett. B345, 558 (1995). 

7. M. Davier and A. Hocker, Phys. Lett. 
B419, 419 (1998). 

8. J.H. Kiihn and M. Steinhauser, Phys. 
Lett. B437, 425 (1998). 

9. J. Erler, Phys. Rev. D59, 054008 
(1999). 

10. S. Groote et al, Phys. Lett. B440, 375 
(1998). 

11. M. Davier and A. Hocker, Phys. Lett. 
B435, 427 (1998). 

12. M. Davier, S. Eidclman, A. Hocker and 
Z. Zhang, Eur. Phys. J. C31, 503 
(2003). 

13. M. Davier, S. Eidelman, A. Hocker and 
Z. Zhang, Eur. Phys. J. C27, 497 
(2003). 

14. K. Hagiwara, A.D. Martin, D. Nomura 
and T. Teubner, Phys. Lett. B557, 69 
(2003). 

15. J.F. de Troconiz and F.J. Yndurain, 
FTUAM-04-02, hep-ph/0402285 (2004). 

16. F. Jegerlehner, DESY-03-189, hep- 
ph/0312372 (2003). 

17. CMD-2 Collaboration 
(R.R. Akhmctshin et al), Phys. Lett. 
B578, 285 (2004). 

18. M. Davier and C. Yuan, Nucl. Phys. 
(Proc. Suppl.) B123, 47 (2003). 

19. V. Cirigliano, G. Ecker and H. Neufeld, 
Phys. Lett. B513, 361 (2001). 

20. V. Cirigliano, G. Ecker and H. Neufeld, 
JHEP 0208, 2 (2002). 



For Publisher's use 



21. Muon g-2 Collaboration (G.W. Bennett 
ef at), Phys. Rev. Lett. 92, 161802 
(2004). 

22. KLOE Collaboration (A. Aloisio et al), 
hep-ex/0407048 (2004). 

23. S. Binner, J.H. Kuehn and K. Melnikov, 
Phys. Lett. B459, 279 (1999); G. Ro- 
drigo, H. Czyz, J.H. Kuhn and M. Szopa, 
Eur. Phys. J. C24, 71 (2002); see 
also the later publications on the subject 
from these authors and collaborators. 

24. BABAR Collaboration (B. Aubcrt et al), 
BABAR-PUB-04/034, hep-ex/0408078 
(2004). 

25. DM2 Collaboration (A. Antonelli et al.), 
Z. Phys. C56, 15 (1992). 

26. V.W. Hughes and T. Kinoshita, Rev. 
Mod. Phys. 71, 133 (1999): T. Ki- 
noshita and M. Nio, hep-ph/0402206 
(2004). 

27. A. Czarnecki and W.J. Marciano, Nucl. 
Phys. (Proc. Suppl.) B76, 245 (1999); 
M. Davier and W.J. Marciano, Annu. 
Rev. Nucl. Part. Sci. 54, 115 (2004). 

28. B. Krause, Phys. Lett. B390, 392 
(1997). 

29. A. Czarnecki, W.J. Marciano and 
A. Vainshtein, Phys. Rev. D67, 073006 
(2003); A. Czarnecki, B. Krause and 
W.J. Marciano, Phys. Rev. Lett. 76, 
3267 (1995); Phys. Rev. D52, 2619 
(1995); R. Jackiw and S. Weinberg, 
Phys. Rev. D5, 2473 (1972); S. Peris, 
M. Perrottet and E. de Rafael, Phys. 
Lett. B355, 523 (1995); M. Knccht et 
al, JHEP 0211, 003 (2002). 

30. K. Melnikov and A. Vainshtein, UMN- 
TH-2224-03, hcp-ph/0312226 (2003). 

31. M. Gourdin and E. de Rafael, Nucl. 
Phys. BIO, 667 (1969); S.J. Brodsky 
and E. de Rafael, Phys. Rev. 168, 1620 
(1968). 

32. CLEO Collaboration (S. Anderson et 
al), Phys.Rev. D61, 112002 (2000). 

33. OPAL Collaboration (K. Ackerstaff et 
al), Eur. Phys. J. C7, 571 (1999). 



34. J. Bijnens and P. Gosdzinsky, Phys. 
Lett. B388, 203 (1996). 

35. S. Ghozzi and F.J. Jegerlehner, Phys. 
Lett. B583, 222 (2004). 

36. M. Davier, hcp-cx/0312065 (2003). 

37. F. Guerrero and A. Pich, Phys. Lett. 
B412, 382 (1997); D. Gomez Gumm, 
A. Pich and J. Portoles, Phys. Rev. 
D62 (2000) 054014. 

38. G. Venanzoni, private communication 
(Aug. 2004). 



